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Abstract 


This article considers a nonlinear inverse problem of the Ostrovsky—Burgers 
equation by using noisy data. Two B-Splines with different levels, the 
quintic B-spline and septic B-spline, are used to study this problem. For 
both B-splines, the stability and convergence analysis are calculated, and 
results show that an excellent estimation of the unknown functions of the 
nonlinear inverse problem. 
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1 Introduction 


In this paper, we consider the Ostrovsky—Burgers equation 
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(ut + bUrae + UUs — AUre)2 = YU+ f(a,t). (1) 


Equation (1) was appeared in modeling internal waves in the ocean or surface 
waves in a shallow channel with an uneven bottom under the effects of the 
interfacial friction (see [11, Chapter 1] and [16, 22]). In (1), the positive 
constants y and a are the rotation and friction coefficients, respectively. The 
function f denotes the external force and 6 is the dispersion coefficient, which 
its sign is related to the type of dispersion. Ignoring the dissipation term uz, 
(1) leads to the Ostrovsky equation 


which was derived by Ostrovsky in 1978 [18] to model weakly nonlinear sur- 
face and internal waves in a rotating ocean; see also [9, 10]. It was also 
demonstrated in [17] that the nonlinear oblique magneto-acoustic waves in 
a rotating plasma can be described by (2). A model of the propagation of 
long internal waves in a deep rotating fluid can be found in [7]. If one con- 
siders the limit of no high-frequency dispersion b = 0, the resulting equation 
is called the Ostrovsky—Hunter equation [4]. When y = 0, (1) turns into the 
Korteweg-de Vries-Burgers equation 


Ut + bUrex + UUg — QUex2 = hace t), (3) 


which is the dissipated version of the KdV equation (a = 0) modeling the 
propagation of weakly nonlinear dispersive long waves in some contexts [19]. 
It is worth noting that despite the similarity of structures of (2) and the KdV 
equation, the Ostrovsky equation, unlike the KdV equation, is nonintegrable 
by the method of the inverse scattering transform [10]. 


In the present paper, we study numerically equation (1) in the domain 
(x,t) € [0,1] x [0,7] with the final time T, the initial condition 


u(x,0) = p(x), xz € [0,1] (4) 
and boundary conditions 
u(0, t) = fi(t), Uz (0, t) = fa(t), te [0, T], (5) 


u(1,¢) = g(t), uz(1,#) = ga(t), Uz (1, t) = g3(t), te (0, T], (6) 


where p(x), gi(t), go(t), g3(t), and f(#,t) are continuous known functions, 
while f(t), fo(t), and the wave amplitude u(z,t) are unknown, which remain 
to be determined. There are different methods to obtain numerical solutions 
of the dissipative Kd V-type equations. 

There are different methods to obtain numerical solutions of the dissipa- 
tive KdV-type equations. Bhatta [2] found numerical solutions of (3) using 
the modified Bernstein polynomials (B-polynomials). In [14], soliton solu- 
tions of the KdVB equation were found by using two numerical methods, 
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finite difference and a semi-analytic method with the Adomian decomposi- 
tion. A spectral method based on Lagrange polynomials was used in [24] 
to approximate solutions of boundary-value problems associated with (3). 
The classical radial basis functions (RBFs) collocation (Kansa) method for 
the numerical solution of (3) was formulated in [12]. A spectral collocation 
method based on differentiated Chebyshev polynomials was elaborated in 
[15] to obtain numerical solutions of (3). Chen and Wu [5] presented the MQ 
quasi-interpolation method to find the numerical solutions of (3). The local 
discontinuous Galerkin method was tested in [23] to study the Kuramoto— 
Sivashinsky equation. 

Here, we consider two numerical methods to find the solutions of (1), 
the collocation method based on septic B-spline basis functions and quintic 
B-spline basis functions. 

It is known that the use of B-splines has many different features and is 
effective in numerical works. One of the most important features is that 
the conditions on the continuity of functions are built-in and have smooth 
interpolation functions. On the other hand, as the support of each B-spline 
is embedded only on a few sub-intervals, the resulting matrix related to the 
discretized equation will be tightly banded. 

Moreover, if one combines with collocation, then the solution procedure 
will be clear and short. In the present work, we provide the combination 
of the finite difference in t and the quintic B-spline collocation method in 
zx, and employ the Tikhonov regularization method to solve the associated 
ill-conditioned system, which gives an efficient explicit solution that has high 
accuracy with the minimal computational effort for the inverse problem (1) 
and (4)—(6). 

This paper is arranged as follows. In Section 2, a description of the septic 
B-splines collocation method, the uniform convergence, and the stability are 
explained. The quintic B-splines collocation method and its convergence and 
the stability are explained in Section 3. The numerical results are presented 
in Section 4, and finally, Section 5 completes this paper with some concluding 
remarks. 


2 Septic B-spline collection method 


To apply the Septic B-spline collection method, region of the solution of 
the problem is restrained over 0 < a < 1. Space interval [0,1] is separated 
into uniformly sized finite elements of length h by the knots z,, like that 
0= 2% <2, <-::< ay =1. Lengths of these finite elements are h = 1-0 — 


N 
(m+1— Lm) form =0,1,2,...,N—1. 
Problem (1)—(6) will be solved with the over-specified conditions 
u(K,t)=hi(t), Us(k,t) =he(t), Uea(k,t) = ha(t), (7) 
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where t € [0,T] and 0 < & < 1 is a fixed number. A spline is a function that 
is piecewise-defined by polynomial functions and possesses a high degree of 
smoothness at the places, where polynomial pieces connect. A B-spline is a 
spline function that has minimal support to the given degree, smoothness, 
and domain partition. We define the septic B-spline B;(x) for 7 = —3(1)N+3 
by the following relations: 


(a — a3~4)", w € [xj-4,%5-2), 
(@= 97-4)’ —8@— 25-2)", & € |wj-3, £j-2), 
(a= 4)' = Sle = 2323) 

4+28(x — x;-2)", x € [Xj;-2, 25-1), 


+28(x x;-2)" 56(a zj;-1)", LE Leas) 
Bj(x) =A" § (wj44— 2)" — 8(aj43 — 2)" (8) 
+28(2j42 — x)" — 56(aj41—2)", 2 € [23,2541), 


(xj44 — 2)" — 8(aj43 — 2)" 


+28(xj;42 — 2)", x E [@j41, Lj 42), 
(aj44 — 2)" — 8(aj43 — 2)", x E [Xj42,2543), 
(ja a)", © € [2543,2544); 
0, otherwise. 


The set 


Q = {B_3(x), B_2(z), B_1(x), Bo(x),..., Bn(2), 
Byii(x), By+o(x), Byis(x)} 


forms a basis for the approximate solution, which will be defined over the 
interval [0,1]. Thus @ = Span(Q) is a subspace of O7[0,1] and w is N + 7- 
dimensional. The values of B; and its derivatives may be tabulated as in 
Table 1. 

Now let U(x,t) € @ be the B-spline approximation to the exact solution 
u(x,t) in the form 


N+3 


U(«,t) = Do qi (t)B,(2), (9) 


j=-3 


where c;(t) are time-dependent parameters determined from boundary con- 
ditions and collocation condition. 


By substituting the trial function (8) into equation (9), the nodal values of 


U, U', U", Us U®), and U“) are obtained in terms of the element parameters 
Cm by 
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Table 1: Values of B;(x) and its derivatives at the nodal points. 


z Uj—-4 | Lj-3 Tj-2 25-1 vj Tj+1 Uj+2 T]4+3 | Vj+4 
B;(«) 0 1 120 1191 2416 1191 120 1 0 
Ba, 0; eg | 2 SP, 0 a 7 Fe 
ay (x) 0 EB us aa =m a ms EB 0 
BY (a) 0 a 160 ee 0 a = see 0 
BM (2) 0 = 0 = 1560 a = 7500 0 = 0 
B®) (x) 0 oe = 10st Pa 0 miei acee a 0 


Um = Cm—3 + 120€m—2 + 1191 em_—1 + 2416¢,, + 1191e, 41 + 1206m+2 + Ems, 


7 
U;, = piv em-8 — 56Cm—2 — 245€m—1 + 2451em41 + 56Cm+2 + Cm+3)s 
Ul = ay 24 15 — 80 15 +24 
m— h2 Cm—3 + Cm—2 + Cm-1 Cm + Cm+1 Cm+2 + Cm+3)s 
Ul = 7; 8 + 19 19 8 
m— ha Cm—3 Cm—2 4 Cm—-1 — Cm4+1 + 8€m+2 + Gineks)> 
840 
ui) = a (em=3 —~ 9Cm—1 +h 16c¢,, a 9Em+41 + Cm+3)s 
2520 
us) — he ( Cm—3 } 4em—2 5Cm—1 t OCm+1 4em+2 t Cm+3)+ 


2.1 Numerical discretization 


Let us consider a uniform mesh (2;,t;) to discretize the region [0,1] x [0,7]. 
More precisely, each (x;,¢;) is defined by x; = th, i = 0,1,2,...,N and 
t; = jk for 7 = 0,1,2,..., where h and k are mesh sizes in the space and 
time directions, respectively (see [8]). Suppose that u” is linearly interpolated 
between two time levels n and n+1 by using the Crank—Nicolson formula and 
then its time derivative is discretized by the forward finite difference formula: 


yrtl tun yrrt 
a ye, (11) 


(et 


where u” = u(z, tn) and wu? = u(x, 0) = p(x). Substituting the above approx- 
imation into equation (1) and discretizing in time variable, we have 


2 
Unt = Yu + OUger — bUerax — Ug, — Uz =F f(a, t), 


sO 
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k 
1 i 1 1 
ugh) — uh = > [ut + yu" + aunts + auitee — (Uties)”*? — (wtes)” 


— (u2)"™* — (u2)” — bur tl, — bu® 


LLLXL LLLELXL 


+f(@tn41) + f(@tn)] 


The nonlinear term is linearized by using the quasi-linearization formula as 
given below: 


Plu ut) = plu ua) + (wrt? — ary OF 4 ant — an QF 


Thus, we have 


ntl = yrtiyn n n+l 


(UUcx) Lu +U Ugg (uter)” 


and ; 
(u2)"* = 2u,"t tu.” — (u2)” : 


Now, we obtain by rearranging the terms that 


Oyu! + bguntt A é3untt a Ogunt} + Osurtt (12) 
= Agu” + ti + O4U re _ Ost ream oT O7( f(a, tnt) ar FG: 


where 


k 7 42 (k 
Q.= = (um, —7), O2=—(1+kut), = 5 (50) 
A 2 

210 (k 840 /k k , 
a= Te (Sa). = Fe (30): a 


By replacing the approximate solution U with u, and using the nodal 
values U and the derivatives of U, we deduce from equation (10) at the knots 
in equation (12) the following difference equation in the variable c: 


Ae Be Bet +4 Creat +p a Bay 4 ees a On, (13) 
= (a, tn) + A(x, tn41); 


where 


A* = 0, — 02+ 03 + 04+ 65, 

B* = 1200, — 5602 + 2403 + 8604, 

C* = 11910, — 24562 + 1503 — 1964 — 96s, 
D* = 24166, — 8063 + 1665, 

E* = 11910, + 24505 + 1563 +1904 — 90s, 
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F* = 1200, +5605 + 2403 — 864, 
G* = 6, + 62 + 63 — 04 + 45, 


(ai, tn) = O6(ci—3 + 120c;_-2 + 1191¢;-1 + 2416c; + 1191ej41 + 120cj+2 + C43) 


7 
() ( ci-3 — 56cj;-2g — 245c;-1 + 2451 ¢;41 + 56cj42 + C443) 


8cj42 + ci+3) 


_ 


210 
as ( \ c—3 — 8cy_-2 +19c;_-1 19¢;4 


840 
95 ( mI ) (4-3 — 9e4—1 + 16c; — 9ei41 + 6443) 


+ O7 f (xi, tn), 


and 
A(2x;,tn41) = 67 f (vi, tn41), 0 < 1 < N. 


System (13) consists of (N + 1) linear equations with (N + 7) unknowns, 


Tr 
(c_3, C_2,C_-1,C0,---,CN,CN+1,€N+4+2; CN+3) 


To have a unique solution to the above system, we are required the over- 
specified conditions (7). Suppose that k = 2,, 1<s< N-—1, then 


u(as,t)=hi(t), us(as,t) = he(t), Uex(as,t) = h(t), (14) 


where t € [0,7]. If we consider m = s in (10), then we have 


hi(tn41) = cRts + 120c7T) + 1191c7*} + 2416cR** + 1191cR tT + 120e7 Ff) + cr ts, 


s s+ 


Te 
ha(tnti) = 5 ch tt — 56cR tS — 245 + 2451cP tT + 56eR > + eR tS), 


AD ps ; ; 
ha(tm+1) = 55 (co Ty + 24eRt> + 15eP tT — 80cR tt + 15ek TT + QdeP 1) + PHS), 
giltnti) = ents + 120cR TI, + 1191chtT + 2416cR*7 + 1191chAT + 120cKH1, + ents, 


7 
g2(tnt1) = 5 (—eQth, — 56cQ TE — 245eK 4 + 2451eN 1 + 56cNt 5 + cK 3) 5 


g3(tn+1) = oa (ently + 24eRtt, + lect — 80cR FT + 15eRt! + 24K + cyt). 
Hence, we derive that 
AC = B, (15) 
is a system of (NV + 7) linear equations with (N + 7) unknowns, where 
A[l,s +1] = A[l,s +7] = A[8,s+1] = A[2,s+7| = A[IN+5,N +7] 
[3,s+7])=A[N+5,N4+1])=A[N+6,N+7| 
[N+7,N+1])=A[N+7,N+7] =—-A[N+6,N +1] =1, 


and 
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N+7,N +2] =A[N +7, N +6] = A[l,s +2] = A[l,s +6] = 120, 
N+7,N +3] = A[IN+7,N +5] = All, s +3] = A[l,s +5] = 1191, 
N+7,N +4] = A[l,s +4] = 2416, 

2,8+6] = —A[2,s +2] = —A[N +6, N +2] = A[N +6, N +6] = 56, 
N+6,N +5] = oe + 6, N +3] = —A[2,5+ 3] = A[2,54 5] = 245, 
3,8 +2) = AB, i 

ee I a 5,N +3) =A[N+5,N +5] =15, 
A[3,s +4] = A[N 


ee 


Therefore, we have 


O ... 1 120 1191 2416 1191 120 1 O 
O ...-1—56—-245 0 245 56 1 O 
O... 1 24 15 -80 15 24 1 O 
A* BX C* D* E* F* G* 0... 
0 A* BY C* D* EX F* G* 0 


A= : Joe, ., a oa om a ’ 
0 50 BY AG? DE: BP oP ae ol) 
0 A’ Be Ot oD" ee ec|* 
0 1 24 15 -80 15 24 1 
0 —1 —56-245 0 245 56 1 
0 1 120 1191 2416 1191120 1 J ining) 
with 


C= oe oo, oft) amee tee — 


(n+1) htt) htt) htt) 
Cy 9 oN41 1 CN 42 7EN43 J > 


T 
a (ay) Bo) Bo Bs Ge tae By Be ees) : 
and 


n ms h i h? 

Bo =hi(tngi), BY = (=) ha(tn4i), BM) = (3) h3(tn41); 
(n) (n) he 

Be’ = A(aijtn) + A(t tnti), Busi = | Zp } galtn+), 


n h n 
B= (7) galtn41), Be, = 71 (tn+41); 
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where 0 <i < N. We notice that the matrix A is ill-conditioned, so we obtain 
the solution of system (15) by using the Tikhonov regularization method. 


2.2 The initial vector c® 


By the initial condition (4) combined with the boundary and over-specified 
conditions (7), the initial vector c° can be written as the following expressions: 


u(xs,to) = co_3 + 120co_, + 1191c2_, + 2416c2 + 1191c8,, + 120c8,. + c2,5 
= hi(to), 


7 
Uz (Ze, to) = ” (—c8_3 — 56c2_ — 245¢2_, + 2451c2,, + 56c2,, + 2) 


= ho(to), 
42 
Une (£5, to) = 53 (co_3 + 24c)_5 + 15c9_, — 80c2 + 15e2,, + 24c2,. + e243) 
= ha(to), 
u(x, to) = ch_3 + 120cP_. + 1191cP_, + 2416c? + 1191c),, + 120c?,. + 2,5 


= P(x), 


where 0 <i< N and 


42 
Ure (tn, to) = - (CN_3 + 24chy_» + 15ch_, — 80), 1 + 15h 41 + 24ch 42 + N43) 


= 9g3(to), 


7 
Uz(@Nn,to) = z (—c)_3 — 56cy_» — 245cR,_, + 2451c0,4, + 56ch 45 + Chi) 


= ga(to), 
u(zn,to) = cy_3 + 120ch_» + 1191ch_, + 2416ch, + 1191ch,) + 120ch 42 + Chas 
= gi (to). 


This yields an (N + 7) x (N + 7)-system of equations of the form of 


AC° = B, (16) 


A[l,s +1] = A[l,s+7| = A[8,s+1] = A[2,s+ 7] = A[IN+5,N +7] 
= AB,s+7]/=A[N+5,N+1]) =A[N+6,N+7| 
=A[N+7,N+4+1)=A[N+7,N +7] =—-A[N+6,N +41] =1, 


and 
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AIN +7,N +2] = AIN +7,N +6] = A[l,s +2] = A[l,s + 6] = 120, 
A[N +7,N +3] =A[N +7, N +5] = Afl,s +3] = All, s +5] = 1191, 
AIN +7,N +4] = A[,s +4] = 2416, 
Al2,s + 6] = —A[2,s +2] = —A[N + 6,N +2] = AN +6, N + 6] = 56, 
A[N +6, N +5] =—A[N +6, N +3] =—A[2, s +3] = A[2, 5 +5] = 245, 
A[3,s + 2] = A[3, s+ 6] = A[N 2] = ALN +5, N + 6] = 24, 
A[3,s +3] = A[3,s+5] = AIN 3] = AIN +5,N +5] =15, 
A[3, s+ 4] = A[N +5, N +4] = —80. 
Thus, we have 
QO... 1 120 1191 2416 1191 120 1 0 
0... —1 —56 —245 0 56 1 «(OO 
QO... 1 24 15 —80 24 1 0 
1 120 1191 2416 1191 120 O sex 
QO 1 120 1191 2416 1191 120 1 0 
A= 
0 1 120 1191 2416 1191 120 1 
0 1 24 —80 15 24 1 
0 —-1 -56 -245 0 245 56 1 
0 1 120 1191 2416 1191 120 1 (N-+7)x(N+7) 
cs hy(to) 
e (7 )ha(to) 
“A (45) ha(to) 
“0 p(xo) 
Os : ‘ udots 
ch p(tw) 
N41 (42)93(to) 
eae (7)92(to) 
CN+3 (N47)x1 91 (to) (N+7)x1 


Finally, the solution of (16) can be obtained by using the Tikhonov regular- 
ization method. 


2.3 Convergence analysis 


Here, we are going to check the convergence of our algorithm. Suppose that 


U(x) = 


N+3 


3 ¢;B;(x) is the B-spline collocation approximation of u(x), 
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where u(x) is the exact solution of (1) with the over-specific conditions (7) 
and initial condition (4). Due to round off errors in computations, we can 
assume that U() is the computed spline for U(x) in such a way that U(x) = 
ake ¢;B;(x), where 

C = (é_3, C_2, C1, Co, cee ,CN; CN4+1) €N+2; én +3). 


Following (15) for U, we have 


where 


B = (hi(tn41), ha(tn41), h3a(tn41), 70, 71,---,7N,93(tn41), 92(tn41), 91(tn41)), 


and 


= O6(ci-3 + 120c¢;-2 + 1191¢;_-1 + 2416c¢; + 1191¢41 + 120¢;42 + Ci+3) 


7 
(7) ( Ci-3 56c;_2 = 245¢;_-1 + 2451 c441 + D6C;+42 + Ci+3) 


210 
64 ( ) ( Ci-3 8c;_2 t 19¢j_1 = 19¢i41 + 8Ci42 + Ci+3) 


840 
— Os (5) (ci_-3 — 9e;-1 + 16¢; — 9¢j41 + Ci+3) 
a 67 (fisted a f (ai, tn41)) : 


Consequently, we have a ” 
A(C —C) =(B- B),. 


The following lemma will be important in our analysis. 


Lemma 1. If {B_3, B_2, B_1, Bo,..., Bn, Bn +1, Bn+2, Bn+3} is the septic 
B-spline, then 


S~> B;(a)| < 7456, x € (0, 1]. (17) 


Proof. At any node x;, we have from Table 1 that 


N+3 N+3 
S> Bj| < S> [Bj] =1+ 120 + 1191 + 2416 + 1191 + 120 + 1 = 5040. 
j=-3 j=-3 


For any point « in each sub-interval [z;_1,2;], we have analogously that 
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N+3 


S> |B; < 7456. 
j=-3 


This proves the lemma. 


To obtain a suitable estimate for the error 
|u(xz) — U(x), 


we need to bound the error terms | w(e) - | and |@ - U(2)| separately. 
To do so, we recall the following theorem. 
Theorem 1. Suppose that g(x) € C®[0, 1] and |g°(x)| < L, for all x € [0,1]. 


Assume that 
A= {0=29 <2, <-:-<ay = 1h}, 


is the partition of [0,1] of step size h. If Sa(a) is the unique spline function 
interpolate g(x) at nodes x,21,...,aN © A, then there exist a constant 
A; <2 such that we have for all a € 0, 1], that 

Ig? (z) — SA@)Iloo SAyLAPY, =f =0,...,7, 


where || - || represents the L°°-norm. 


Proof. See [1] for a proof. 


Now first we find a bound on ||B — Blloo- Following [20, Theorem 7] and 
applying Theorem 1, there exists W = Wz > 0 such that 
|e-Bl| <w (ju@ = O(«)| nS \U'@) = 0’ (2)| Es lu") ~ 6" (x) 
+|U" (2) - 8" (@)| + ]U @) - @))) 
<WLAgh® + WLAh" + WLAgh® + WLA3h® + WLgh*. 


Thus we can get . 
|B -— Blloo < Wih*, (18) 


where 
W, = WL(Agh* + Ah? + Agh? + Agh + Aa). 


Since the matrix A is ill-conditioned, we have from the Tikhonov regulariza- 
tion solution that 


(c-@) =[4r4+a(x)"] 4" (BB), 
after taking the D°-norm, 
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a T = a 
Jo-2]_s|[ata+a (a) "| atl pal) < wan 
_ (19) 
where 
. =i 
W2=W, [ATA +a (x) re Ay 
Now we observe that 
= N+3 
U(x) —U(«) = So (ej — &)Bi(2). 
1=-3 
By taking the L°-norm and using (17) and (19), we have 
a N+3 N+3 
|u@-F@]|_ =| Me@-aHB@] <ie-al, | Oo BO) 
i=—3 oe i=-—3 
< 7456W2h’4. 
(20) 


Theorem 2. The time discretization process (11) that we use to discretize 
equation (1) in time, is convergent of the first order. 


Proof. See [21] 


Theorem 3. Let u € C®[0, 1] be an exact solution of (1) such that 


oun 


sel <x, 


for all x and t. If U(a,t) is the numerical approximation by our method of 
u, then 


||u(x) — U(2)llo < O(k +h’). 
Proof. We have from Theorem 1 that 
||u(a) — U(x) loo < Ap Lh. (21) 


Equations (20) and (21) imply that 


u(x) - U(a)Il,, < |[ule)-G(@) |) + || F@)-Ve)|| 
< AoLh® + 7456W2h* = Thi, 


Co 


where T = \oLh* + 7456W2. Therefore we deduce from Theorem 2 that 
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llu(x) — U(e)llo S o(k +h’), 


where p is a finite positive constant independent of k and h. This completes 
the proof. 


2.4 Stability analysis 


In this section, we will investigate the stability by applying the Von-Neuman 
stability analysis. Without loss of generality and for the sake of simplicity, 
we assume in equation (1) that f(z,y) = 0. On the other hand, as we 
have linearized the nonlinear term uu, by considering u as a constant k, in 
equation (1), then the equation can be rewritten by 


(ut + bUrae + kiUe — QUcx)2 = YU. (22) 


Substituting the approximation (11) into equation (22), by discretizing in 
time we have that 


—Pyyu®* + uth! + Putt — Pautt) + Putt 


LLL LLLXL 


= Pyu” + ur, ~ Pour, ay Py ~~ 4) 
If E =T, then 
42 210 840 
Pp=7T7, Po= p27 P3 = ye Th a ye TO: (23) 


In terms of unknown time parameters c;, the equation can be written from 
(10) by 


n+1 n+1 n+1 n+1 n+1 n+1 n+1 
L1C; "3 + Lg; "9 + X3C; "7 + Lac; + U5) + LEC o + L7Ci 3 


—_ n n n n nt+1 n+l n+1 
= £€gC;_3 + LoCji_9g + Vi0c;_ 1 + V11¢C7 + C16) + Lisees + Lia; 
+ + +3 


where 


ay =—-P,-1+P,4+ P3+ Pa, 

is, = S1D0P) — 56 +24) 4 BPs, 

v3 = —1191P, — 245+ 15P) — 19P3 — 9P,, 

a4 = —2416P, — 80P2 + 1664, (24) 
ie 1 101 Pp ae 4 RPS EOP: = GP, 

ve = —120P, + 56 + 24P, — 8P3, 

v7 =—P,+1+ Po— P3+ Fa, 
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tg=P,-1-P,-—P3-—F,, 
tg = 120P, — 56 — 24P, — 8P3, 
tio = 1191 P, — 245 — 15P,+19P3 + 9P4, 
111 = 2416P, + 80P, — 1664, 
@12 = 1191 P, + 245 — 15 Py — 19P3 4+ 9P4, 
%13 = 120P, + 56 — 24P, + 8P3, 
ti4= Pip +1—Po+ Ps — Py. 


Now we substitute cf”, = AC” exp(im@h), where 1 = /—1, A is amplitude, h 
is the step length and £ is the mode number. We get 


Xy4i% 


Sas er 


where 


X1 = (vg + ®14) cos(38h) + (xg + 213) cos(28h) + (x19 + 212) cos(8h) + xe, 
X22 = (%7 + 41) cos(3BhA) + (x2 + x6) cos(28h) + (x3 + #5) cos(Bh) + va, 

Y, = (414 — &g) sin(38h) + (#13 — xg) sin(28h) + (a12 — £19) sin(Bh), 

Y2 = (a7 — #1) sin(38h) + (ve — x2) sin(28h) + (v5 — x3) sin(Bh). 


(25) 
Hence, we have from (25) and (24) that X, = X9. Setting 
p = 2sin(3Gh) + 112 sin(26h) + 490 sin(3h), (26) 
yp = P3(2sin(3Gh) + 16 sin(28h) — 38 sin(Bh)), 
we have that 
Yi=pt+y, Yo=p-y. 
To get the stability, we need to prove |¢| < 1, that is, 
XP + (e+)? < XZ +(p—- 9). (27) 


To this end, it is enough to show that 4py < 0. If we rewrite (26) by 


p = 8sin(Gh)(cos?(Bh) + 28 cos(Bh) + 61), 
yp = 8P3sin(8h)(cos(bh) + 5)(cos(Gh) — 1), 


then we obtain 
App = 256P3 sin? (8h)(cos? (8h) + 28 cos(Bh) + 61)(cos(Bh) + 5)(cos(8h) — 1). 


By the positivity of the constants, P3 > 0. The functions sin?(8h), cos? (Bh)+ 
28 cos(3h) + 61, and cos(Gh) + 5 are positive, but cos(Gh) — 1 is negative, 
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so 4pp < 0. Hence, we obtain that the seventh order B-spline collocation 
method is unconditionally stable. 


3 Quintic B-splines collocation method 


In this section, we are going to solve the inverse problem (1) by a new mod- 
ification of the quintic B-splines collocation method with the over-specified 
conditions 


WN 


u(K,t) = hi(t), 
Uz (kK, t) = ho(t), 48) 


where t € [0,t7] and 0 <& <1 is a fixed number. 


We define the quintic B-spline B;(x) for 7 = —2(1)N +2 by the following 
relation: 


(ewes), x € [rj_3,£j_2), 
(f—Gya)? — Ole — 2-5)”, © € |[xj~2, 25-1), 
(x — 23-3)” — 6(@ — x32)? + 15(x — 2j-1)°, 2 € [2j7-1, 24), 
B(a) =h~? & (aj43 — 2)? — 6(aj42 — 2)? + 15(x541 — 2), x € [25,25 41), 
(343 — 2) — 6(r;42 — 2)”, BE [Para Bes), 
(2343 — 2)”, x € [@j42,%j43), 
0, otherwise. 
(29) 


We note that the set 
Q = {B_2(zx), B_i(zx), Bo(z),..., Bn (x), By 4i(z), Bn+2(z)}, 


forms a basis for an approximate solution that will be defined over the 
interval [0,1]. Hence, the set @ = Span((Q) is a subspace of O?{0,1] and w 
is N + 5-dimensional. The values of B; and its derivatives are tabulated in 
Table 2. Let U(a,t) € @ be the B-spline approximation to the exact solution 


u(x,t) in the form 
N+2 


U(a,t) = S° cj (t)B;(e), (30) 


j=-2 


where c;(t) are time-dependent parameters determined by the boundary and 
collocation conditions. 


Substituting the trial function (29) into (30), the nodal values of U, U’, 
U",U, and U™) are obtained in terms of the element parameters Cm by 
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Table 2: Values of B;(x) and its derivatives at the nodal points. 
w Uj-3 | Vj-2 | Vj-1 | Vj | Vj+1 | Vj+2 | V+3 
B(x) 0 1 26 66 26 1 0 
7 5 50 —50 5 
Bi (a) 0 h 7 OO ae |) 0 
Wa 20 10 —120 40 20 
By (2) 0 h2 h2 h2 h2 he 0 
60 —120 120 —60 
By(a) | 0 ne he 0 | ae 0 
(4) 120 —480 720 —480 120 
B (x) 0 hé h4 he he ht 0 


Um = Cm—2 + 26Cm—1 + 666m + 26Cm+1 + Cm+2; 
5 


UL = 7 (-em—2 = 10Gm—1 + 10 en + Oma), 
Ul = a 2 6 2 
mn = 75 (Cm—2 + 2m—1 — 6m + 2em41 + Cm+2), (31) 
yi _ ea 9) 2 
m= res —Cm—2 + 2Cm—1 — 2Cm4i + Cents 
120 
uM) = sa (em-2 — Aem—1 + 6Cm — 4Cm+1 + Cm+2)- 


3.1 Numerical discretization 


To discretize the equation, similar to the previous section, we use (31) in 


2 
Unt = Yu + QUree bUrrex Wer Ua f(a, ae 
and obtain 


k 
1 1 1 1 
Up s —Up = 5 [yume + yu" + aupge + Upe, — (Uuer)"** — (utter)” 


2 
— aay 7 (uz)" a nee - DU si 


+f (x, tn41) + f(z, tn)] 3 


The nonlinear term is also linearized by rearranging terms as given below: 


burt + Oguntt + O3untt = ic ssailea + Osun) = (32) 
Ogu” or Up + O4U roe = Oat ree + O7( f(x, tn41) a5 f(z,tn)), 


where 
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: =: : 20 (/k 60 (/k 
1 = 5 (Ure — 7) O2 = 7 (1+ kuz), oa = 3 (Su) a= 9 (50), 


Again, by replacing the approximate solution U with u, and using the 
nodal values U and the derivatives of U, we obtain from equation (31) at the 
knots in equation (32) the following difference equation in the variable c: 


Atchty + Bch! + Otcht) + D*chitl + Bt chy = H(2i,tn) + Hai, tn41), 
(33) 

where 

A* = 6, — 62 + 03 + 04 + 85, 

B* = 260, — 1062 + 2603 — 204 — 46s, 

C* = 660) — 603 + 665, 

D* = 2606, + 202 + 203 + 204 — 40s, 

EX = 0, + 09 + 63 — 64+ 95, 


A(a;, bin) = 6 (ci_-2 + 26c;-1 + 66¢; + 26c;41 + Ci+2) 


5 
a (7) (—ci-2 = 10¢;_1 + 10cj41 + Ci+2) 


60 
+ 04 @ (—ci-2 + 2ci-1 — 2ci41 + Ci+2) 
120 
Os ( 74 ) (ci_-2 Ac;-1 + 6c; — 4cj41 + Ci+2) 
+ O7f (ri; tn), 


and 
A(a;,tn41) => 07 f (@i,tn41), 0 < 1 < N. 


System (33) contains (IV + 1) linear equations with (NV + 5) unknowns, 


T 
(c_2, C_1,C0,---,CN,CN+1, CN+2) 
To have a unique solution to the above system, we are required the over- 
specified condition (28). Assume that « = x,, 1 < s < N—1. Equation (6) 
holds and moreover 


u(as,t) =hi(t), Uux(as,t) = ho(t), 


where t € [0,7]. If we consider m = s in (31), then we have 
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Ai(tn41) = ets + 26cRt} + 66c241 + Wert! + erty, 
5 
ho(tn41) = i (erty + 1l0cRt} — Lock t} — etf5), 


(tr41) 

gi (tn41) = CNa + 26K + 66cRT! + 26M + chit, 
( )= = (cht, + 10cKT, — l0cK TY + ee) - 

Consequently, 

AC =F; (34) 


is a system of (N +5) linear equations with (N +5) unknown functions where 


Afl,s+1] = A[l,s+5] = A[2,s+1)=A[N+4,N +45] = A[2,54+5] 
= AIN+4,N +1) = AIN+5,N +5] =A[N +5, ee 
AIN +5,N +2] = A[N +5,N +4] = A[l,s +2] = A[l,s +4] = 
AIN +5,N +3] = A[l,s +3] = 66, 
A 
A 


N+4,N +2] = A[2,s+2] =-A[N+4,N +4] =—A[2,s +4] =10, 
2,6+3)=A[N+4,N +3] =0. 


Thus, we can write the matrix A as 


O ... 1 26 66 26 1 O 
0... 1 10 0 -10 1 O 
A* B* C* D* EX 0... ... 
0 A* B* C* D* E* 0 


A= ; 
0 A* B* C* D* E* 
0 1 10 O -10 1 
0 Ts 26 86.926. Joe ois 
with 
n+1 n+1 n+1 n n+1 n+1 n+1 
C= as her of I sitege: aera ee Some deny" ; 


de 
BBB e BT sccghl eB on ey a 


n n h n 
BS =M(tny), BO = (Z) ha(ingi), BS”) = H(ai,tn) + A (ai, try); 


n h nm 
Bee (Z) galtn41), Bo. = 7 (tn+41); 
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where 0 <i < N. We note that A is an ill-conditioned matrix, so we obtain 
the solution of system (34) by using the Tikhonov regularization method. 


3.2 The initial vector c® 


From the initial condition (4) the boundary and over-specified conditions (28) 
the initial vector c° can be obtained as the following expressions: 


u(xs, to) = Co_5 + 26c2_, + 66c2 + 26c2,, + C249, 
Uz (25, bp = > («2 a) an We) 1 1024 1 Ou) = aa ha(to), 


) 
) 
u(xj, to) = co» + 26c)_, + 66c? + 26c?,, +c... = p(zi), 0<i<N, 
) 
) 


= h > (ae a + 10ch_1 — 10ch 41 + Cyi2) = g2(to), 
Cr 


og + 26c4,_ 1 + 66c%, + 26c%4, + N42 = gi(to). 
This gives an (N + 5) x (N + 5)-system of equations of the form of 
= B, (35) 


where 


> 


1,s+1)=All,s+5] = A[2,s+1) = A[N+4,N +5] = A[2,5+5], 
=A[N+4,N4+1)=A[N+5,N +5] = A[N+5,N4+1]=1, 


AIN +5,N +2] =A[N +5,N +4] = A[l,s +2] = A[l,s + 4] = 26, 
A[N +5,N +3] = A[l,s +3] = 66, 
AIN +4,N +2] = A[2,s +2] =—A[N +4,N +4] =—A[2,s +4] = 10, 


> 


2,s+3] = A[IN+4,N +43] =0, 
Thus, we have 


O... 1 2666 26 1 O 
0... 1 10 0 —10-1 0 


126 66 26 1 O ...... 
0 1 26 66 26 1 O 


1 26 66 26 1 


0 1 10 0 —10 1 
0 6a 1 26 66 26 1 


(N+5) x (N+5) 
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cae hi(to) 
oy (#)ha(to) 
ay p(xo) 
Oe 2 p= 
cy p(rn) 
CN+1 ()g2(to) 
CnN+2 (N+5)x1 gi(to) (N+5)x1 


Finally the solution of (35) can be obtained by using the Tikhonov regular- 
ization method. 
3.3 Convergence analysis 


Similar to the Convergence of the previous part, first we need to recall the 
following lemma. 


Lemma 2. The B-splines {B_2, B_1, Bo,..., Bn, Bn +1, Bn+2} satisfy the 
following inequality: 


S~ B;(x)| < 186, (0<a2<1). (36) 


Proof. We have from Table 2 at any node 2;, that 


N+2 N+2 
Y Bix >. (B= 0. 
i=—2 i=-—2 


Similarly, we have for any point x in each sub-interval [x;_1,x;] that 


N+2 


S> |B;| < 186. 


j=—2 


Theorem 4 ([6, 13]). Suppose that f(x) € C®[0,1] and |f®(x)| < L, for all 
x € [0,1]. Assume that 


A={0=29 <4 <-:-<ay = 1h}, 
is the partition of [0,1] of step size h. If Sa(x) is the unique spline function 


interpolates f(x) at nodes x,71,...,2~ € A, then there exist a constant , 
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such that we have for all x € [0,1] that 


|F@) 7 S4(2)|| SL, gO Noid. 
where || - || represents the L°-norm. 


By using the ideas of [20] and Theorem 4, there exists W = W; > 0 such 
that 


uy 


|B = B| <W ((U@) = G(2)| 4: lu’@) = 0" (2)| a: lu"(@) ~ 8" (2) 
+|0" (2) - 8" (@)| + |U@) - 6 ~@))) 
< WLAgh® + WLAh® + WLAgh* + WLAgh? + WLAgh?. 


Thus we can get 
|B- Bl < Wih?, (37) 


where W, = WL (Aoh* + Ah? + Agh? + Ash + Ya) - Since the matrix A is 
an ill-conditioned matrix, we have from the Tikhonov regularization solution 
that 


(0-0) =[474+0(n)" 2] 4" (B-B). 


After taking the infinity-norm we find from (37) that 


_fi T (2)\7 sl T 
Ig él|_< [4 A+a(R ) re A 


|B - 3 <Wyh?, (38) 


where 
W2=W, [ATA + a(R)? R2)-1 AT | 


co 


Now we see that 


N+2 


U(x) —U(«) = S~ (ci, — &) Bila). 


1=—2 


Using (36) and (38), we have after taking the L°-norm that 


= N+2 N+2 
U@ -F@_ =| @-wB@) <Me-al. | XO Be) 
i=—2 ag i=-2 
< 186W2h?, 


(39) 
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Theorem 5. Let u(x,t) € C®[0, 1] be the exact solution of (1) such that 


6 
Assume that U(a,t) is the numerical approximation by our methods, then 
||u(z) — U(x) loo S$ O(k + h?). 
Proof. We have from Theorem 4 that 
||u(x) — U(2) loo < AoLh®, (40) 


Equations (39) and (40) imply that 


lulz) — U(a)Ilo < |[ule)-G(@) |] + | F@)-ve)|| 
< ApLh® + 186W4h? = Th, 


co 


where TY = \oLh* + 186W2 and U(z,t) is the approximate solution by our 
present method. Then by Theorem 2 and the above considerations, we have 


||u(x) — U(z)loo S$ p(k +h”), 


where p is a finite positive constant independent of k and h. This completes 
the proof. 


3.4 Stability analysis 


Similar to the previous section, we use the Von-Neuman stability analysis. 
We assume again in equation (1) that f(x,y) = 0. On the other hand, as we 
have linearized the nonlinear term uu, by considering u as a constant k, in 
the equation (1), then the equation can be rewritten by 


If we substitute the approximation (11) into the equation (41) and discretize 
it in time, we obtain that 
—Pyu"*! + qr + Poutt? _ Pzunt} + Pyutt 
= Pu" + uh — Pout, + Paut, — Pau 


n 
LLLXL* 


If £ =T, then we have 
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60 20 120 
Pi =TY;, P3 = hee Py, _ 727k Py = ya Te (42) 


In terms of unknown time parameters c;, the equation can be written from 
(31) by 


* n+1 * n+1 * n+l * n+1 * n+l 
os he ee Ree a es 
_ pen eon en * n+l * n+l 
= F'ce g + Ghee + Bee + Pe + Ics 


where 


A= —P,—-14+2)—P, +2, 

B* = —26P, + 10+ 2P) + 2P3 —4pu, 

C* = -66P, — 6P, + 6Pr, 

D* = =96P, — 10 +9; —2P; =4P,, 

B= Py 14 Fe Pp Py, (43) 
F*=P,+1—-P,+P3—Pi, 

G* = 26P, +10 —2P, —2P;+4P,, 

H* = 66P, + 6P) — 6Fi, 

I* = 26P, —10—2P) +2P3 +4Pu, 


J* =P, —1—P,— P3— Py. 


If we substitute c?, = AC” exp(imGh), where h is step length, A is amplitude, 
i = /—1, and £ is mode number, then we get that 


X1,+1Y, 
, = N+ y’ 
where 
X, = (£* + J*)cos(28h) + (G* + I*) cos(8h) + A", 
Xq = (A* + E*)cos(26h) + (B* + D*) cos(8h) + C, 
Y, = (J* — F*) sin(2Gh) + (I* — G*) sin(Gh), 
Y2 = (E* — A*)sin(26h) + (D* — B*) sin(Bh). 


We have from (43) that X2 = —X, and 


Y, = sin(28h)(—2 — 2Ps3) 4 
Y2 = sin(28h)(—2 + 2Ps3) 4 


+ sin(Bh)(—20 + 4P3), 
- sin(Bh)(—20 — 4P3). 


To show the stability, we need to prove |¢| < 1, that is, 


X24 V7 < X?4+Y¥? 
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We prove that Y? — Y? < 0, but we see that 
Y? — Y? = 256P3 sin?(8h)(cos(Bh) — 1)(cos(Bh) + 5). 


As all arbitrary constants are positive, then P; > 0. Moreover, the functions 
sin?(8h) and cos(Gh) + 5 are positive but cos(8h) — 1 is negative, so Y? — 
Y? <0. Hence, we obtain that the quintic B-spline collocation method is 
unconditionally stable. 


4 Numerical result and illustrations 


In this section, we will give some examples to investigate the applicability of 
our methods that we described in the previous sections. In general, there are 
two error sources for an inverse problem. The first one is the unavoidable bias 
deviation, known as deterministic error, and the second one is the variance 
caused by the enhancement of error estimation, known as stochastic error. 
The global impact of deterministic and stochastic errors is studied in the 
mean squared error or the total error as in [3]. Next, by considering the total 
error S defined by 


NIB 


N 
Sf, = ro Sen) = fad? , 


s=1 


NI 


N 
Sf. = a SUF) ~ flea)? , 


s=1 
we will compare the exact and the approximate solutions, where f; and fo 
are exact values, ff and f3 are the estimated values, and the number of 
estimated values is denoted by N. 
We should remark that there are no previous numerical results related to 
the nonlinear inverse problem (1) and (4)—(6) in the literature to compare our 
results, and one may use the other numerical methods to study this problem. 


Example 1. In our first example, we consider the nonlinear inverse problem 
(1) and (4)-(6), where a = 1, b= 5, and y = 3 with the initial data 


u(z,0) = sin(x), 
and the external force 
3 
f(a, t) =cos(t +a) + 5 sin(t + x) + cos*(t + 2) — sin?(t +2). 


An exact solutions of this problem is u(a,t) = sin(a + t) with 
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u(0,t) = filt) = sin(t), 
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uz(0,t) = fo(t) = cos(t). 


Table 3: Comparison between exact and numerical solutions of Example 1 for u(0.1, t), 
|u(0.1,t) — u*(0.1, t)|, by the Quintic and Septic B-spline methods with N = 30, 50, 100. 


N=30 N=50 N=100 

t Quintic Septic Quintic Septic Quintic Septic 
0.1 4.8e — 04 4.2e — 06 8.7e — 04 6.5e — 06 1.4e — 03 8.5e — 06 
0.2 7.3e — 04 4.4e — 06 9.0e — 04 4.9e — 06 1.8e — 03 3.7e — 08 
0.3 1.0 — 03 1.6e — 06 1.7e — 03 1.9e — 06 2.3e — 03 2.0e — 06 
0.4 1.4e — 03 1.5e — 05 2.2e — 03 6.9e — 06 3.0e — 03 5.5e — 06 
0.5 1.7e — 03 1.0e — 05 2.7e — 03 1.8e — 08 3.5e — 03 5.3e — 08 
0.6 1.9e — 03 1.le — 05 2.9e — 03 3.2e — 06 3.8e — 03 3.3e — 06 
0.7 2.3e — 03 1.8e — 05 3.3e — 03 4.le — 05 4.3e — 03 2.2e — 06 
0.8 2.6€ — 03 3.0e — 06 3.7e — 03 2.0e — 06 4.7e — 03 3.6e — 06 
0.9 2.8e — 03 2.4e — 06 3.9e — 03 3.6e — 06 4.9e — 03 1.4e — 07 

1 3.0e — 03 1.3e — 06 4.2e — 03 2.0e — 06 5.2e — 03 2.1e — 06 
Su | 5.9758e — 05 | 1.0042e — 06 | 8.7071le — 05 | 1.0093e — 06 | 1.131le — 04 | 1.0143e — 06 


Table 4: 


N = 30, 50, 100. 


Comparison between exact and numerical solutions of Example 1 for 
fit), |fi() — ff (|, by employing the Quintic and Sep 


ic B-spline methods with 


N=30 N=50 N=100 

t Quintic Septic Quintic Septic Quintic Septic 
0.1 0.9e — 04 1.9e — 05 1.le — 04 5.9e — 05 1.2e — 04 1.9e — 06 
0.2 2.4e — 04 2.2e — 05 2.4e — 04 3.5e — 06 2.8e — 04 7.3e — 08 
0.3 3.9e — 04 3.le — 05 3.4e — 04 8.5e — 08 4.6e — 04 5.0e — 08 
0.4 6.4e — 04 3.4e — 05 5.le — 04 9.4e — 06 6.6e — 04 1.4e — 07 
0.5 3.2e — 04 2.8e — 05 6.2e — 04 1.3e — 05 4.0e — 04 2.0e — 07 
0.6 1.0e — 03 2.7e — 05 7.2e — 04 3.6e — 06 4.6e — 04 1.8e — 07 
0.7 1.2e — 03 2.3e — 05 8.7e — 04 4.7e — 06 5.2e — 04 4.5e — 08 
0.8 1.3e — 03 2.2e — 05 9.6e — 04 2.4e — 07 5.6e — 04 3.3e — 08 
0.9 1.4e — 03 1.4e — 05 1.0e — 03 3.7e — 06 6.le — 04 4.3e — 08 

1 1.6e — 03 1.5e — 06 l.le — 03 3.7e — 06 6.5e — 04 4.6e — 07 
Sy, | 3.1833e — 05 | 1.1633e — 06 | 2.3106e — 05 | 1.1575e — 06 | 1.338le — 05 | 9.8944e — 07 


Tables 3-5 show the total error S' for some values of N for each method. 
In Figures 1 and 2, we plotted the total error S as a function of N. Also 
we took « = N~!, T = 1, At = 0.001, and the noisy data (input data + 


0.001 x rand(1)). 
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Table 5: Comparison between exact and numerical solutions of Example 1 for f(t), 
fo(t) — f3(t)|, by the Quintic and Septic B-spline methods with N = 30, 50, 100. 
N=30 N=50 N=100 
t Quintic Septic Quintic Septic Quintic Septic 
0.1 4.9e — 03 1.9e — 03 3.6e — 04 6.3e — 03 1.2e — 02 7.9e — 06 
0.2 1.2e — 02 2.3e — 03 8.le — 03 4.9e — 03 1.9e — 02 5.2e — 04 
0.3 1.7e — 02 3.5e — 03 1.3e — 02 3.7e — 03 2.6e — 02 4.2e — 04 
0.4 2.3e — 02 4.1le — 03 2.2e — 02 2.5e — 03 3.3e — 02 1.0e — 04 
0.5 2.8e — 02 3.7e — 03 2.8e — 02 2.4e — 03 4.0e — 02 1.5e — 04 
0.6 3.4e — 02 4.0e — 03 3.4e — 02 9.le — 04 4.6e — 02 8.3e — 05 
0.7 3.9e — 02 3.6e — 03 4.le — 02 1.2e — 03 5.2e — 02 4.6e — 04 
0.8 4.3e — 02 3.9e — 03 4.6e — 02 3.7e — 04 5.6e — 02 1.7e — 05 
0.9 4.7e — 02 3.3e — 03 5.le — 02 3.le — 04 6.le — 02 5.le — 05 
1 5.2e — 02 3.5e — 03 5.5e — 02 3.0e — 04 6.5e — 02 2.3e — 05 
S;, | 9.6783e — 04 | 7.7698e —05 | 0.0012 8.392le—05 | 0.0013 1.5842e — 06 


Plot of variation error of u when N=30 Plot of variation error of u when N=50 


Plot of variation error of u when N=100 


Figure 1: Plots of variation u(a,t) of Example 1 by using the Quintic B-spline method 
when N = 30,50, 100. 


Example 2. In this example, we consider the inverse problem (1) and (4)- 
(6), where a = 1, b=3 and y = 0.5. The initial data is 
u(a,0) = sech(z), 


while 
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Plot of variation error of u when N=30 Plot of variation error of u when N=50 


er) 


Plot of variation error of u when N=100 


0s 0.8 


02 


Figure 2: Plots of variation u(x,t) of Example 1 by using the septic B-spline method 
when N = 30,50, 100. 


oe 142sinh(t— x)? sinh(t — x)? — 6sinh(¢ — x)? + 2sinh(¢ — x)? 
eS ~~ cosh(t — 2)3 cosh(t — «)4 
120sinh(t — x)* 23. Bsinh(t— x) — 1 


cosh(t — x) cosh(t— 2x) '  cosh(t—2)? ’ 


is the external force. An exact solution of this problem is u(,t) = sech(a—t) 
such that 


u(0,t) = f(t) = sech(—t), 


sinh(t) 

Uz (0, t) = f(t) _ cosh(t?) | 

Tables 6 and 7 give the numerical results for conditions u(0,t) and u,(0,t), 
respectively. To explain the accuracy of our methods, some illustrations are 
presented in Figures 3-6. Table 8 shows the numerical values of u(x,t) at 
point z = 0.1. Figures 7 and 8 also present the comparison between the 
numerical and exact values of u(a,t). Also we take « = 0.03,T = 1, At = 
0.001 and the noisy data (input data + 0.001 x rand(1)). 
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Table 7: Comparison between exact and numerical solutions of Example 2 for f(t), 
fo(t) — fZ(t)|, by the Quintic and Septic B-spline methods with N = 100. 


Quintic B-spline Septic B-spline 

: RO | FO | AO-HOL|L_ AO] fO_| AO-HOI 
0.1 0.099172 | 0.021154 | 7.801727e — 02 | 0.099172 | 0.102857 | 3.685425e — 03 
0.2 0.193493 | 0.118409 | 7.508375e — 02 | 0.193493 | 0.197758 | 4.265252e — 03 
0.3 0.278678 | 0.208712 | 6.99660le — 02 | 0.278678 | 0.282586 | 3.907888e — 03 
0.4 0.351456 | 0.292363 | 5.909310e — 02 | 0.351456 | 0.355666 | 4.210337e — 03 
0.5 0.409814 | 0.356937 | 5.287678e — 02 | 0.409814 | 0.413920 | 4.105471e — 03 
0.6 0.453029 | 0.414228 | 3.880057e — 02 | 0.453029 | 0.457373 | 4.344355e — 03 
0.7 0.481503 | 0.454315 | 2.71877le — 02 | 0.481503 | 0.486087 | 4.583750e — 03 
0.8 0.496500 | 0.479423 | 1.707765e — 02 | 0.496500 | 0.500676 | 4.175531le — 03 
0.9 0.499829 | 0.495907 | 3.922220e — 03 | 0.499829 | 0.500025 | 1.958806e — 04 

1 0.493554 | 0.499845 | 6.290468e — 03 | 0.493554 | 0.492724 | 8.303327e — 04 
Sf, - - 0.0018 - - 1.3344e — 04 

Execution time (second) 40.079962 45.367252 
Condition Number of Matrix A 1.4332e + 20 4.1816e + 19 


Table 6: Comparison between exact and numerical solutions of Example 2 for f1(t), 


fit) — ff (|, by the Quintic and Septic B-spline methods with N = 100. 
Quintic B-spline Septic B-spline 
é A®_ | FO | A®-FO| AO] AO_] AO-AO 
0.1 0.995021 -997416 | 2.394818e — 03 | 0.995021 | 0.995007 | 1.369526e — 05 
0.2 0.980328 | 0.982549 | 2.220868e — 03 | 0.980328 | 0.980312 | 1.635728e — 05 
0.3 0.956628 -958612 | 1.984447e — 03 | 0.956628 | 0.956613 | 1.456806e — 05 
0.4 0.925007 | 0.926788 | 1.780935e — 03 | 0.925007 | 0.924991 | 1.596319e — 05 
0.5 0.886819 -888307 | 1.488184e — 03 | 0.886819 | 0.886804 | 1.534899e — 05 
0.6 0.843551 | 0.844734 | 1.183566e — 03 | 0.843551 | 0.843534 | 1.641813e — 05 
0.7 0.796705 | 0.797546 | 8.406108e — 04 | 0.796705 | 0.796696 | 9.782520e — 06 
0.8 0.747700 -748193 | 4.927829e — 04 | 0.747700 | 0.747694 | 6.334825e — 06 
0.9 0.697795 | 0.698057 | 2.627784e — 04 | 0.697795 | 0.697795 | 3.326556e — 08 
1 0.648054 -647956 | 9.857496e — 05 | 0.648054 | 0.648057 | 2.488186e — 06 
Sp, - - 5.0435e — 05 - - 5.0706e — 07 
Execution time(second) 40.079962 45.367252 
Condition Number of Matrix A 1.4332e + 20 4.1816e + 19 


Figure 3: Comparison between exact and numerical values of (a) u(0,t) and (b) 
ux(0,t) of Example 2 by the Quintic B-spline method. 
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Figure 4: Comparison between exact and numerical values of (a) u(0,t) and (b) ua (0, t) 
of Example 2 by the septic B-spline method. 


Plot of variation u 


a asurface u(x,t) 
# Ulx,t) 


0.5 


Figure 5: Comparison between exact and numerical values of u(x,t) of Example 2 by 
the quintic B-spline method. 


Plot of variation u 


surface u(x,t) 
+ _U(xst) 


Figure 6: Comparison between exact and numerical values of u(x,t) of Example 2 by 
the septic B-spline method. 
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Table 8: Comparison between exact and numerical values of u(0.1,t) of Example 2 
by the Quintic and Septic B-spline methods with N = 100. 


Quintic B-spline Septic B-spline 

t u(0.1,t) | u*(0.1,t) | Ju(0.1,¢) — u*(0.1,¢)] | u(0.1,t) | u*(0.1,t) | Ju(0.1,¢) — u*(0.1, t)| 
0.1 | 1.000000 | 0.996042 3.958350e — 03 1.000000 | 1.000008 8.493153e — 06 
0.2 | 0.995021 | 0.991131 3.889610e — 03 995021 | 0.994971 4.942389e — 05 
0.3 | 0.980328 | 0.976786 3.542139e — 03 980328 | 0.980355 2.729571e — 05 
0.4 | 0.956628 | 0.953011 3.616572e — 03 956628 | 0.956594 3.413866e — 05 
0.5 | 0.925007 | 0.922351 2.656163e — 03 925007 | 0.925047 3.991245e — 05 
0.6 | 0.886819 | 0.884390 2.429210e — 03 .886819 | 0.886839 1.991189e — 05 
0.7 | 0.843551 | 0.841742 1.808390e — 03 843551 | 0.843521 2.948822e — 05 
0.8 | 0.796705 | 0.795713 9.927146e — 04 796705 | 0.796733 2.729211e — 05 
0.9 | 0.747700 | 0.747000 7.002802e — 04 -T47700 | 0.747675 2.541404e — 05 

1 | 0.697795 | 0.697668 1.267156e — 04 .697795 | 0.697777 1.739320e — 05 
Su - - 9.9114e — 05 - - 8.2254e — 07 


Figure 7: Plot of variation u(0.1,t) of Example 2 by the Quintic B-spline method with 
N = 100. 


“exact 


Figure 8: Plot of variation u(0.1,t) of Example 2 by the septic B-spline method with 
N = 100. 


Example 3. In our last example, we solve the inverse problem (4) and (4)- 
(6), where a = 1, b= 2, y = 0.5 and the initial data is given by 
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u(x, 0) = cos(z), 


The force function is 


5 cos(t — 
= cost ®) sti 


f(a,t) x) — cos(t — 2)? + sin(t — x)’, 

The function u(x, t) = cos(x — t), is an exact solution of this problem. Thus 
u(0,t) = fi(t) = cos(—t) and u(1,t) = fo(t) = —sin(-t). Tables 9 and 10 
give the numerical results for u(0,t) and u,(0,t), respectively. To explain the 
accuracy of our methods, some illustrations are presented in Figures 9-12. 
Table 11 shows the numerical values of u(x,t) at point « = 0.1. Figures 13 
and 14 also present the comparison between the numerical and exact values 
of solution u(z,t). Also we take « = 0.03, T = 1, At = 0.001, and the noisy 


data (input data+0.001 x rand(1)). 


Table 9: 
applying 


Comparison between exact and numerical values for f1(t) of Example 3 by 
the Quintic and Septic B-spline methods with N = 100. 
Quintic B-spline Septic B-spline 


t+ 


fi(é) 


fi) 


lA) — AO! 


fi) 


Fi) 


lA) -— Ol 


oO 


0.995004 


0.996818 


1.813881le — 03 


0.995004 


0.994960 


4.446296e — 05 


0.980067 


0.981816 


1.749680e — 03 


0.980067 


0.980097 


3.086339e — 05 


= 


0.955336 


0.957039 


1.702169e — 03 


0.955336 


0.955366 


2.917672e — 05 


0.921061 


0.922675 


1.613838e — 03 


0.921061 


0.921150 


8.917805e — 05 


0.877583 


0.879143 


1.56091 le — 03 


0.877583 


0.877252 


3.307205 — 04 


So 


0.825336 


0.826760 


1.424698e — 03 


0.825336 


0.824996 


3.394589e — 04 


0.764842 


0.766191 


1.349309e — 03 


0.764842 


0.764432 


4.104127e — 04 


So 


0.696707 


0.697931 


1.224282e — 03 


0.696707 


0.696210 


4.966064e — 04 


>] G0) Na] @] er] AY] bo] tf 


0.621610 


0.622693 


1.082930e — 03 


0.621610 


0.621069 


5.406920e — 04 


= 


0.540302 


0.541248 


9.459846e — 04 


0.539825 


4.770495e — 04 


Shr 


4.6373e — 05 


0.540302 


3.8840e — 06 


Execution time(second) 


40.562310 


41.0305 


30 


Condition Number of Matrix A 


1.4332e + 20 


1.3648e + 19 


Table 10: Comparison between exact and numerical values for f2(t) of Example 3 by 
applying the Quintic and Septic B-spline methods with N = 100. 
Quintic B-spline Septic B-spline 


t+ 


fo(t) 


f(t) 


lfo(t) — fa) 


fo(t) 


f3 (4) 


[fa(t) — fa) 


0.099833 


0.038846 


6.09871le — 02 


0.099833 


0.103929 


4.095972e — 03 


oO 


0.198669 


0.139510 


5.915962e — 02 


0.198669 


0.194482 


4.186958e — 03 


0.295520 


0.238522 


5.699784e — 02 


0.295520 


0.292183 


3.337169e — 03 


0.389418 


0.335087 


5.433175¢e — 02 


0.389418 


0.378289 


1.112973e — 02 


0.479426 


0.427911 


5.151492e — 02 


0.479426 


0.466016 


1.340928e — 02 


0.564642 


0.516943 


4.769963e — 02 


0.564642 


0.546380 


1.826272e — 02 


9S 


0.644218 


0.599703 


4.451437e — 02 


0.644218 


0.692126 


4.790832e — 02 


0.717356 


0.677664 


3.969196e — 02 


0.717356 


0.774241 


5.688534e — 02 


o 
&] G0] Na] @] Or] GA] Go] tof 


= 


0.783327 


0.748561 


3.476550e — 02 


0.783327 


0.846131 


6.280393e — 02 


0.810553 


3.091790e — 02 


6.059395e — 02 


0.841471 


0.0016 


0.841471 


0.902065 


5.1323e — 04 


Execution time (second) 


40.562310 


41.0305 


30 


Condition Number of Matrix A 


1.4332e 4 


+ 20 


1.3648¢ 4 


+ 19 
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Figure 9: | Comparison between exact and numerical values of (a) u(0,t) and (b) 


uz(0,t) of Example 3 by the Quintic B-spline method. 


Figure 10: Comparison between exact and numerical values of (a) u(0,t) and (b) 
uxz(0,t) of Example 3 by the septic B-spline method. 
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Figure 11: Comparison between exact and numerical values of u(x,t) of Example 3 
by the Quintic B-spline method. 
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by the septic B-spline method. 


Comparison between exact and numerical values of u(x,t) of Example 3 


Table 11: Comparison between exact and numerical values of u(0.1,t) of Example 3 
by the Quintic and Septic B-spline methods with N = 100. 


Quintic B-spline Septic B-spline 

t u(0.1,t) | u*(0.1,t) | |u(0.1,¢) — u*(0.1,¢)| | u(0.1,t) | u*(0.1,t) | |u(0.1,¢) — u*(0.1, t)| 
0.1 | 1.000000 | 0.996544 3.455663e — 03 1.000000 | 0.999902 9.781401le — 05 
0.2 | 0.995004 | 0.991605 3.399369e — 03 995004 | 0.995019 1.462144e — 05 
0.3 | 0.980067 | 0.976735 3.331920e — 03 980067 | 0.980016 5.030518e — 05 
0.4 | 0.955336 | 0.952079 3.257457e — 03 955336 | 0.955342 5.215045e — 06 
0.5 | 0.921061 | 0.917956 3.104844e — 03 921061 | 0.921182 1.214016e — 04 
0.6 | 0.877583 | 0.874601 2.981736e — 03 877583 | 0.877510 7.290696e — 05 
0.7 | 0.825336 | 0.822554 2.781252e — 03 825336 | 0.825410 7.484661e — 05 
0.8 | 0.764842 | 0.762264 2.577809 — 03 764842 | 0.764753 8.910966e — 05 
0.9 | 0.696707 | 0.694359 2.347544e — 03 .696707 | 0.696628 7.884439e — 05 

1 | 0.621610 | 0.619528 2.082224e — 03 621610 | 0.621652 4.249717e — 05 
Su - - 9.5839e — 05 - - 2.8097e — 06 


Figure 13: Plot of variation u(0.1,t) in Example 3 by the Quintic B-spline method 


with N = 100. 
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Figure 14: Plot of variation u(0.1,t) in Example 3 by the Septic B-spline method 
with N = 100. 


5 Conclusion 


In this paper, we have successfully employed the septic B-spline and the 
quintic B-spline method to estimate unknown boundary conditions in an 
inverse problem related to the Ostrovsky—Burgers equation (1) and (4)-(6). 
We have discussed the convergence rate of our methods and have shown 
the rate of the septic B-spline method is O(k + h*), while O(k + h?) is 
the convergence rate for the quintic B-spline method. The stability of both 
methods was investigated. By comparing the numerical results, we showed 
that the accuracy and stability of the septic B-spline method are more than 
the ones for the quintic B-spline method. Since the associated coefficient 
matrix in the septic B-spline method and quintic B-spline method are usually 
ill-conditioned, we have used the Tikhonov regularization method to obtain 
a stable numerical approximation of the solution. The results are collected 
by using the MATLAB 7.10(R2016a), tested on a personal computer with 
Intel(R) Core(TM)2 Duo-CPU and 4-GB RAM. 
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